Crossover from 0(3) to 0(4) behavior in weakly frustrated 

ant if er r omagnet s . 



Andrey V. Chubukov 1,2 and Oleg A. Starykh 3 'Q 
1 Department of Physics, University of Wisconsin, Madison, WI 53706 
2 P.L. Kapitza Institute for Physical Problems, Moscow, Russia 
3 University of California at Davis, Davis, CA 95616 

(today) 

Abstract 

We consider an anisotropic version of the CP 1 model which describes frus- 
trated quantum antiferromagnets with incommensurate spin correlations. We 
extend the two-component spinon field, describing lattice spins, to the M- 
component complex vector, and show, in the 1/M expansion, that for arbi- 
trary small incommensurability longitudinal and transverse stiffnesses tend 
to the same value as the system approaches the quantum critical point. For 
physical spins (M = 2), this yields 0(4) critical behavior. However, if the 
spin structure is commensurate, the longitudinal stiffness is identically zero. 
In this case, the critical behavior is the same as in 0(3) sigma model. We show 
how the critical exponents interpolate between 0(3) and 0(4) values near the 
transition. We also show that the competition between these two fixed points 

leads to a confinement-deconfinement transition at a finite temperature. 
PACS: 67.50-b, 67.70+n, 67.50Dg 
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In this communication, we address a specific issue concerning the nature of zero- 
temperature quantum phase transitions in quantum antiferromagnets: whether there is a 
qualitative change in the critical behavior when one adds a frustration to Heisenberg antifer- 
romagnet and makes spin configuration non-collinear. We argue that the critical behavior 
of collinear and non-collinear antiferromagnets is qualitatively different, even when local 
spin configuration differs from a Neel state by an arbitrary small amount. As a quantita- 
tive measure of this difference, we consider the ratio of the fully renormalized transverse 
and longitudinal stiffnesses, p±/p\\. In a collinear antiferromagnet, the longitudinal spin 
stiffness is identically zero because the Neel ordering is described by just one vector of anti- 
ferromagnetism, n, and a rotation of n around its equilibrium direction is not a symmetry 
transformation. On the other hand, in the ordered non-collinear state, one needs two unit 
vectors to describe the equilibrium spin configuration. In this case, a rotation of one vector 
around another is a legitimate symmetry transformation which is broken in the ordered 
state. As a result, there are two finite stiffnesses: one for two equivalent transverse spin- 
wave modes, and another for the longitudinal spin-wave mode. The action for non-collinear 
antiferromagnets can be written either in terms of the SO (3) rotation matrix, or in terms 
of two-component complex spinon field. Azaria et al. performed one-loop RG studies of the 
5*0(3) action [||] and found that as the system approaches the critical point, the ratio of the 
fully renormalized stiffnesses tens to unity. Sachdev, Senthil and one of us used the spinon 
description and have studied the critical behavior of frustrated antiferromagnets under the 
assumption that the bare transverse and longitudinal stiffnesses, p\ and pjj, are close to each 
other [0. They also found that as the system approaches the critical point, the ratio of the 
fully renormalized stiffnesses tends to unity. At criticality, p± = pn, and the symmetry of 
the underlying action is enlarged from SU(2) x U(l) to 0(4) (see below). In this paper, 
we extend the spinon approach to an arbitrary ratio of the bare stiffnesses. We will show 
that the 0(4) critical behavior holds for arbitrary pjj/ p° x < 1. We cannot say at the moment 
about how far the 0(4) behavior extends to a region where pjj S> p\. 

Our point of departure is the partition function for a frustrated antiferromagnet written 
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in terms of spinon fields, Z = J Dz*Dzexp [-S], where 
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Here z is a two-component complex spinor field subject to constraint z* a z a — 1, /i — r, x, y, 
and 7 = (pjj — p°±)/p°±- For simplicity we choose units where h — 1 and set both spin-wave 
velocities to unity. 

This effective action can be explicitly derived from semiclassical microscopic considera- 
tions P,|J and the general macroscopic approach of Ref f|. Note however that z quanta 
are not Schwinger bosons. The relation between z and the underlying spins is more com- 
plex and involves incommensurate ordering momentum @||]. The action in ([!]) is invariant 
under global SU(2) spin rotation, and is also invariant under certain type of lattice trans- 
formations. As shown in 0], this lattice symmetry is in essence identical to a lattice U(l) 
symmetry. We will thus refer to the total global symmetry of the effective action in (|l|) as 
SU(2) x U(l). It is essential however, that the U(l) gauge symmetry z(r, r) — > z(f, r)e l ^ r ' T ^ 
is broken provided pu ^ 0, i.e., 7 > — 1. It is only present at 7 = — 1 in which case the action 
in ([!]) describes collinear antiferromagnets. In this latter case, the description in terms of z 
quanta is equivalent to the description in terms of Schwinger bosons: the same action as in 
(H) is obtained if one introduces the U(l) gauge invariant Schwinger boson decomposition 
n a = into the partition function of the 0(3) n— field model 7< 1 limit of the 

action (P was studied in detail in Ref. 0. Here we focus on a region near 7 = — 1. 

To perform 1/M expansion, we need to generalize the action to large M. We generalize 
the doublet z to the M-component complex vector, rescale the z field to z — > z/yM (such 
that z* a z a = 1, a = 1,2...M), and introduce the coupling constant g = M/2p\. We further 
introduce the Hubbard-Stratonovich vector gauge field to decouple the quartic term, and 
introduce a constraint into the action using the integral representation of the 5— function. 
We then obtain 




(2) 
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where we introduced r = — (1 + 

The most straightforward way to compute the ratio of the fully renormalized stiffnesses, 
which we will follow, is to perform calculations in the ordered state at T = 0. This state 
is realized for g smaller than the critical coupling g c . Assume that the first component of z 
is condensed. We then write z = (a,ir a ), and represent a as a sum of the condensed part, 
VMuo, and fluctuation a around it, a = \^Ma + a. It is also convenient to introduce pairs 
of real variables instead of complex variables a and n a : a = x + i 1 *}, = 02a- 1 + i<p2a, 
and rescale gauge field as — > (1 + r)A^. Substituting these expressions into the action 
we find 



(d, X ) 2 + (<V?) 2 + £ (dM 2 - 2VMa oV d,A, 
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This decomposition of the z field implies that the variables <p a describe 2M — 2 transverse 
fluctuations, r\ is a variable for a longitudinal mode, and x describes fluctuations in the 
direction of the condensate. 

Our first goal is to integrate out fluctuations of x, 77, and (p a and to obtain the effective 
action for collective variables A and A. The integration over the longitudinal and transverse 
fluctuations yields contributions to the effective action which are linear in M, whereas the x 
field contributes only a subleading, 0(1), term which can be safely neglected in the leading 
order calculations in 1/M. Performing the Gaussian integration over fluctuating fields and 
using the constraint equation at M = 00, we obtain the effective action for the gauge 
fluctuations in the form Sa = (M/2) j d 3 q U^ u (q) A^(q)A u (—q), where 

d 3 k 1 f d 3 k (2k fl + q^)(2k u + q u ) 
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U ^ q) = 2 J(2^¥-JT2^ 



(2nfk 2 J (2tt) 3 k 2 (k + q) 2 

+ 25^±±^± - MS* . (4) 
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The first two terms in (^) are the components of the polarization operator of the 0(3) sigma 
model, n^(g), which should be massless due to the gauge invariance of the latter. Using 
the Pauli-Willars regularization, we obtain II^,(g) = (8^ — ^f 1 )^- Collecting all terms in 
(||) and inverting the result we find for the gauge field propagator 
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Notice that the longitudinal part of the propagator appears only due to incommensurability 
||. This is a direct consequence of the fact that incommensurability breaks the gauge 
symmetry of the action @. 

The propagator of the constraint field, IT _1 (g), can be calculated in a similar way. We 
found that, to leading order in 1/M, 11(g) is independent on r and has the same form as in 
§: n- 1 (g) = 8g 2 /(g + 16a 2 /^). 

We proceed now with the calculations of the stiffnesses at M = oo. The two stiffnesses 
can be extracted from the long-distance behavior of the propagators of the transverse and 
longitudinal fields: G~7^(q) = p±q 2 , G~ l (q) = P\\q 2 - The computation of these propagators 
is straightforward. Transverse fields <fi a do not directly couple to the condensate, and all 
corrections to the free-particle propagator have relative 1/M smallness. Hence full M = oo 
propagator coincides with the bare one G,p a (q) = g/(2q 2 ), i.e., p± = 2/g. The longitudinal 
(rj- field) propagator is however different from G^ a already at M = oo because of the coupling 
term oc ^/Maor/d^A^ in (|3]). This term leads to the finite M = oo correction, 

G,h) = ^ + Mai q ~fD,M = ^ (6) 

so that p|| = g^+o-^ • We emphasize that Eq.(^) is an exact M = oo result, not an expansion 
around the free-particle expression. 

For the ratio of the fully renormalized stiffnesses we then obtain 



This is the key result of the M = oo consideration. We see that as long as the longitudinal 
stiffness is finite (i.e., 7 > —1), the ratio of the fully renormalized stiffnesses approaches as 
the system moves to the critical point, o"o ~~ > 0. At criticality, 7 = 0, and the renormalized 
action ([]]) reduces to that for the 0(2M) sigma model. At the same time, if pj| = (i.e., 
7 = —1), then 7 — 7 = — 1, and the renormalized action retains the symmetry of the 
isotropic CP M ~ l model. Alternatively stated, for pn = 0, which is the case for a collinear 
antiferromagnet, the system is in the basin of attraction of the CP M ~ X fixed point. However, 
an arbitrary small amount of non-collinearity drives the system to the 0(2M) fixed point. 

Our next goal is to compute 1/M corrections to the ratio of stiffnesses. These corrections 
come from self-energy diagrams which involve exchange by fluctuations of both the constraint 
and the gauge field. The computational steps are rather involved, but conceptually are 
similar to those discussed in Refs. |7|Q]. For r 3> (Tq we obtained 

a? = ~ fi-d^i) («) 
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where L = log(g c — g)/g c , and o"o is related to g c — g as in the 0(2M) model 0, a\ = 
(1 — g/gc) 1 _4 /( 7I " 2a/ ). This is merely a consequence of the fact that the propagator of the 
constraint field does not depend on r. For the ratio of the stiffnesses we then obtain 7 = 
(7/(1 + 7)) (1 — g /g c ) 1+16 ^ 3n M \ For 7<C 1, this reproduces the result of Ref. @. We see 
that the 1/M corrections only speed up the flow to the 0(2M) fixed point [fLl |. For r < Uq, 
the expression for 7 is rather involved and we refrain from presenting it. 

For completeness, we also computed critical exponents in the 1/M expansion. For 7 = — 1 
(i.e., r = 0) spinon correlation function and order parameter possesses r\ = — -fL, v = 
1 + -?Tf, 2/5 = 1 — -2m- The actual spin susceptibility is a convolution of two spinon 
fields, and it has different critical exponents fj, u, etc. We computed the spin susceptibility 
to order 1/M and found 77 = 1 — ^f|j, u — v, (3=1 + • The result for fj has been 
reported by us previously [0. 

For 7 > — 1 (i.e., r > 0), the gauge field acquires a mass, and the self-energy terms 



associated with the exchange of gauge field fluctuations are no longer singular. We have 
checked that the spinon fields now possess 0(2M) exponents J7[: r\ = * M , v — \ — 
3 ^2 M , 2/5=1 — For spin-spin correlation function and magnetization (at arbitrary 

7) we reproduced the results of Ref ||: (3 = 1 + O (j^z) , v — 1 — 3 ^ M , fj — 1+ ■ 

So far we were discussing the zero-temperature critical properties of incommensurate 
antiferromagnets. As the exponent fj is finite for both 7 = — 1 and 7 > — 1, the staggered 
static spin susceptibility, x(q)> a ^ criticality possesses a branch-cut singularity independent 
on whether the system is in the basin of attraction of 0(3) or 0(4) fixed points. From 
this perspective, the difference between 0(3) and 0(4) critical behavior at T = is only 
quantitative but not qualitative one. In general, however, the descriptions of the system in 
terms of n— field and in terms of spinons are fundamentally different: in the first case the 
excitations necessary possess integer spin while in the latter one can have excitations with 
either integer or half-integer spin depending on whether spinons are confined or not. We now 
show that the presence of two fixed points at T = gives rise to a confinement- deconfinement 
transition within the disordered region. Consider for definiteness the renormalized - classical 
(RC) regime, g < g c ||. Here the ground state is ordered, but an arbitrary small temperature 
leads to the restoration of the spin rotational symmetry which implies that one can no 
longer distinguish between transverse and longitudinal fluctuations. Consider first collinear 
antiferromagnet, 7 = —1. Then the critical behavior is governed by the 0(3) exponents. As 
the isotropic CP 1 model is isomorphic to the 0(3) sigma- model, one should obtain the same 
result within the n-field and the spinon description. We recently demonstrated that this 
is indeed the case. Despite the fact that effective action (HD yields a branch cut behavior 
of x(q) a ^ the mean-field (M — > 00) level, the gauge field fluctuations, which appear at the 
1/M level, confine spinons into pairs with integer spin. These bound states of spinons yield 
poles in x(<z)> an d the long-distance behavior of spin correlators totally consistent with the 
fi-field description. The branch-cut to pole transformation is a direct consequence of the 
gauge-invariance of the action (JI|) at 7 = — 1. Unbroken gauge invariance leads to a gap less 
gauge field fluctuations, which give rise to the unbounded long-range confining potential 
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between spinons. 

Let us consider now what happens at 7 > — 1. We have shown in that the static 
staggered susceptibility x(o) is proportional to the Green's function ty(q,x = 0) of the 
effective inhomogeneous Schrodinger equation 

, d 2 



dx 2 



+ V(x)+5 2 )V(x)=5(x) } (9) 



where 5 2 = ^- + Mq, and tuq is the gap in the spinon spectrum. At M = 00, the confining 
potential V(x) vanishes, and x{q) nas a branch-cut singularity at 5 = 0. At finite M, 
V(x) is given by the regularized Fourier transform of the transverse part of the gauge field 
propagator || 

m T f°° dk^ „_ Wi _ ikx . 



V{X) = ~W loo 2^ Dtrans {k){1 ~ 6 ] - (10) 



In the RC regime, we found that in the small momentum limit, D trans {k) is 

(11) 



Dtrans(k) ~ ^ 



k 2 T 2r n_ ' 



+ — 



127rmQ ' g 

The appearance of the mass term 2r/g in D trans (k) makes the would-be-confining po- 
tential between spinons short-ranged: V(x) = (Q^m^/ Mm^) (1 — e~^ mA ), where = 



moy2AiTr/gT. Our key observation is that the gauge field mass, m^, is inversely propor- 
tional to the temperature. At high temperatures rriA is small, and the potential is linear 
in \x\ upto large scales leading to a strong binding of spinons ||. However, as T — > 
rriA — > 00, and V(x) reduces to a constant in which case no bound states exist. There 
exists therefore a critical value of the gauge field mass when the attraction between spinons 
becomes too weak to bound them into pairs. Careful analysis of the homogeneous ver- 
sion of Eq.(|9]) performed by Campostrini and Rossi [[J shows that bound states disappear 
when (6vr/M) 1/3 (m /m A ) < C, where C = 0(1) is a numerical constant. Confinement- 
deconfinement temperature is then 

T* = 8C 2 (6n/M) 1/3 p\r. (12) 
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Above T* there is a confinement, and the staggered static spin susceptibility x(o) nas a P°l e 
singularity, Xconf(q) = 2+ m 2 {A — > as T — > T*), which translates into x(r) oc r _1 / 2 e~ rm at 
large distances. Below T*, spinons are deconfined, and x(<z) has only a branch-cut singularity, 
Xdeconf(q) ~ (<? 2 + t^ 2 )" 1 ^ 2 which implies that at large distances x( r ) °t r~ x e~ Tm . There is 
therefore a real change in the behavior of the physical observable at T*, however, we do not 
expect that there will be any changes in the thermodinamic quantities at this temperature. 

We also caution that our low-energy analysis is restricted to the neighborhood of the CP 1 
fixed point, where the gauge field mass is small at T ~ T*. The fate of the transition line 
at larger r is unknown simply becase one cannot use the Hubbard- Stratonovich decoupling 
when the gauge field mass becomes comparable to the upper cut-off of the theory Jl2| . 



The results of the present paper are in agreement with Ref || and with the results of 
other authors. Transformation from free spinons at r 7^ to the confined ones at r = 



was observed by Wiegmann [|10] in the exact solution of the 2D classical 0(3) problem. 
Deconfinement of spinons in the T = quantum-disordered phase with incommensurate spin 
correlations, due to the appearance of the gauge field mass, was discussed by Sachdev and 



Read JT3|. Classical ci-dimensional versions of the actions ([I]) and (|3|) were recently studied by 
Azaria et al. |T2j by the renormalization group and 1/M analysis in all dimensions between 2 
and 4. One of their key findings is that all models with 7 7^ — 1 are asymptotically equivalent 
at long distances to the 0(2M) model. This is in complete agreement with our results. 
They also argued that though the low-energy physics is controlled by 0(2M) fixed point, 
the high-energy (short distance) behavior for weak incommensurability is still controlled by 
the the CP 1 fixed point - this is also consistent with our result that confinement persists at 
high temperatures. They, however, did not explicitly discuss the confinement-deconfinement 
transition for the spin susceptibility. 

To conclude, in this paper we analyzed the crossover from the CP 1 = 0(3) to the 
0(4) critical behavior in a model of weakly frustrated quantum antiferro magnet. We have 
shown, in the 1/M expansion, that for an arbitrary small bare longitudinal stiffness, the 
system flows away from the QP M ~ X fixed point towards the 0(2M) fixed point. The same 
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type of critical behavior was found previously for the case where the bare longitudinal and 
transverse stiffnesses were close to each other ||. It is therefore likely that the 0(2M) 
critical behavior (i.e., the 0(4) behavior for the physical case of M — 2) holds at least for all 
< Pj_- An unresolved issue is whether the 0(4) behavior holds for an arbitrary large ratio 
of pjj/pj-5 or t nere is a crossover to a different kind of critical behavior with possible binding 
of spinons. This last issue is interesting on its own grounds but it is also possibly related 
to the controversy surrounding the critical behavior of stacked triangular antiferromagnets 
which, as numerical studies indicate, possess critical exponents different from the 0(4) ones 



14j| . The ratio of the stiffnesses in 2D triangular antiferromagnets is not known exactly: 



noninteracting spin- wave calculations yield p?J Pj_ = 2 but first 1/S corrections are not small 
and substantially reduce this ratio ]15| . 



We also have shown that that the competition between the two zero-temperature fixed 
points leads to a confinement- deconfinement transition for static spin susceptibility at a 
finite temperature. Numerical studies of this phenomena are very desirable. 
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